XP-00206693U io..oS-\^ r 

Celt. Vol. 90. 797-507. August 22, 1997. Copyright 31997 by Colt Press ^Jtj - $0 ^. ■ - 

Ancient Missense Mutations in a New Member ~ 
of the RoRet Gene Family Are Likely to Cause 
Familial Mediterranean Fever 



The International FMF Consortium* 



Summary 

Familial Mediterranean fever (FMF) is a recessively 
inherited disorder characterized by dramatic episodes 
of fever and serosal inflammation. This report de- 
scribes the cloning of the gene likely to cause FMF 
from a 115-kb candidate interval on chromosome I6p. 
Three different missense mutations were identified in 
affected individuals, but not in normals. Haplotype and 
mutational analyses disclosed ancestral relationships 
among carrier chromosomes in populations that have 
been separated for centuries. The novel gene encodes 
a 3.7-kb transcript that is almost exclusively expressed 
in granulocytes. The predicted protein, pyrin, is a 
member of a family of nuclear factors homologous to 
the Ro52 autoantigen. The cloning of the FMF gene 
promises to shed light on the regulation of acute in- 
flammatory responses. 



Introduction 

The hereditary periodic fever syndromes are a relatively 
recently recognized group of inherited inflammatory dis- 
orders that are characterized by episodic, self-limited 
bouts of fever accompanied by unexplained arthritis, 
sterile peritonitis, pleurisy, and/or skin rash (reviewed 
in Kastner, 1996). To date, three clinical diseases have 
been identified in this class: the recessively inherited 
familial Mediterranean fever (FMF, MIM249100) and 
Dutch type periodic fever (MIM260920, also known as 
the hyperimmunoglobulinemia D syndrome), and the 
dominantly inherited familial Hibernian fever (MIM1 42680). 
In addition, there are a number of other related but 
clinically distinguishable conditions (Gertz et al., 1987; 
Karenko et al M 1992) that have been described in case 
reports but have not yet been organized into specific 
nosologic entities. 

Of the three well -defined hereditary periodic fever 
syndromes, FMF is by far the most common and has 
attracted the most investigative interest. FMF is seen 



'The International FMF Consortium is 
comprised of the following groups: 
Group 1: 

Ivona Aksentijevich.'t Michael Centola.'t 
Zuoming Deng, 1 * Raman Sood/t 
James E. Balow, Jr., 1 Geryl Wood,' 
Nurit Zaks, 1 Elizabeth Mansfield/ 
Xiang Chen,' Shlomit Eisenberg/ 
Anil Vedula, 1 Neta Shafran, 1 Nina Raben,' 
Elon Pras/ Mordechai Pras/ 
and Daniel L. Kastner't 
'Arthritis and Rheumatism Branch 
National Institute of Arthritis and 

Musculoskeletal and Skin Diseases 
Bethesda, Maryland 20892-1820 
department of Medicine C 
HHeller Institute of Medical Research 
Chaim Sheba Medical Center 
Tei-Hashomer 
Israel 52621 
Group 2: 

Trevor Blake, 4 Andreas D. Baxevanis/ 
Christiane Robbins/ David Krizman/ 
Francis S. Collins/ and P. Paul Liu 4 
'Laboratory of Gene Transfer 
^Genome Technology Branch 
'Laboratory of Cancer Genetics 
National Human Genome Research Institute 
Bethesda, Maryland 20892 
Group 3: 

Xiaoguang Chen/ Mordechai Shohat,* 
Melanie Hamon/ Tanaz Kahan/ 
Andrea Cercek/ Jerome I. Rotter/ 
and Nathan Fischet-Ghodsian 7 



'Departments of Pediatrics and Medical 

Genetics 
Cedars-Sinai Medical Center 
Los Angeles, California 90048-0750 
^Department of Medical Genetics 
Beilinson Medical Center 
Petah Tiqva 
Israel 49100 
Group 4: 

Neil Richards, David A. Shelton, 

and Deborah Gumucio 

Department of Anatomy and Cell Biology 

University of Michigan 

Ann Arbor, Michigan 48109-0616 

Group 5: 

Yuji Yokoyama. Marie Mangelsdorf, 
Annette Orsborn, and Robert I. Richards 
Department of Cytogenetics and Molecular 
Genetics 

Adelaide Women's and Children's Hospital 
North Adelaide 
South Australia 5006 
Group 6: 

Darrell O. Ricke, Judith M. Buckingham, 

Robert K. Moyzis, Larry L. Deaven, 

and Norman A. Doggett 

Center for Human Genome Studies 

Los Alamos National Laboratory 

Los Alamos, New Mexico 87545 



tThase authors contributed equally to this work and 

are listed in alphabetical order. 

*ro whom correspondence should be addressed. 



Cell 
798 



primarily in individuals of non-Ashkenazi Jewish, Arme- 
nian, Arab, and Turkish background. In these popula- 
tions, the carrier frequency has been estimated to be 
as high as 1 :5 (Rogers et al. t 1989; Daniels et al., 1995; 
Yuval et al., 1995), raising speculation that heterozy- 
gotes may have some selective advantage, perhaps 
manifested by increased resistance to an as yet uniden; 
tified infectious agent. Depending on poorly defined he- 
reditary and environmental factors (Pras et al., 1982), 
homozygotes often develop progressive systemic amy- 
loidosis from the deposition of the acute phase reactant 
serum amyloid A (SAA). Prior to prophylactic treatment 
with colchicine, FMF amyloidosis was a significant 
cause of renal failure and death in the Middle East. In 
the absence of an accurate diagnostic test, FMF patients 
in countries where the disease is less prevalent often 
experience years of attacks and several exploratory sur- 
geries before the correct diagnosis is made and ade- 
quate treatment with colchicine is begun, 

FMF attacks are characterized by a massive influx of 
polymorphonuclear leukocytes (PMNs) into the affected 
anatomic compartment, inciting factors are unclear, 
and, even off treatment, patients can sometimes inexpli- 
cably go for long periods without symptoms, only to 
relapse into a frequent pattern of attacks at some later 
time. At the biochemical level, patients have variously 
been reported to have abnormal levels of a C 5a inhibitor 
(Matzner and Brzezinski, 1984), neutrophtl-stimuiatory 
dihydroxy fatty acids (Aisen et al., 1985}, and dopamine 
p-hydroxylase (Barakat et al., 1988), but none of these 
findings has fed to the identification of the genetic basis 
of FMF. 

In the absence of other functional clues, we undertook 
a positional cloning approach to identify the FMF gene 
(designated MEFV). Initial linkage studies placed MEFV 
on chromosome 1 6p (Pras et aL, 1 992) in a panel of non- 
Ashkenazi Jewish families. Linkage to this region has 
subsequently been confirmed and extended to all four 
major affected ethnic groups {Shohat et al., 1992; Pras 
et al., 1994; French FMF Consortium, 1996), although 
differences ^n severity among subpopulattons has sug- 
gested the possibility of allelic heterogeneity (Pras et al., 
1997). Early studies demonstrated a substantial founder 
effect in the Moroccan Jewish population (Aksentijevich 
et at., 1 993), and as markers closer to MEFV have been 
identified, it has become apparent that a high percent- 
age of all North African Jewish FMF carrier chromo- 
somes (93% in our hands) are descended from a com- 
mon ancestor {French FMF Consortium, 1 996; Levy et 
al., 1 996; Balow et al., 1 997). This same ancestral haplo- 
type was observed in one third of the Armenian carrier 
chromosomes we have studied and in about 15% of 
Iraqi Jewish FMF chromosomes. 

Using publicly available polymorphic markers, we nar- 
rowed the candidate region to an approximately 1 Mb 
interval lying between the polycystic kidney disease 
{PKD1) and tuberous sclerosis (TSC2) genes on the teto- 
meric end, and the CREB-binding protein (CREBBP) 
gene on the centromeric end (Sood et al„ 1996; see 
Figure 1). The extensive physical maps constructed 
around these genes (Dackowski et al M 1996; Giles et al., 
1997) did not, however, extend into the MEFV region, 
and large insert YACs mapping to this interval (Doggett 



et al., 1995) were found to be unstable. We therefore 
established a highly redundant contig comprised of 
small-insert YACs, P1s, PACs, BACs, and cosmids that 
completely spans the interval (Sood et al., 1997). Using 
an independent panel of non-Ashkenazi Jewish families 
and newly developed microsatellites from Genethon, the 
French FMF Consortium (1996) identified markers -600 
kb apart defining the MEFV interval by intrafamilial re- 
combinants. By comparing FMF carrier chromosomes 
in their panel with the ancestral haplotype noted above, 
the French Consortium found "historical recombinants" 
that further narrowed the candidate interval to a —250 
kb YAC. Our own data defined a 285 kb candidate inter- 
val by observed intrafamilial recombinants and an —200 
kb minimal region by historical recombinants (Balow et 
al., 1997). Both of these sets of data are consistent in 
their placement of the FMF gene. 

We now report the cloning of the FMF gene from 
this 200 kb interval. Our strategy involved developing a 
detailed transcript map from this region by a combina- 
tion of exon amplification, direct cDN A selection, single- 
pass sequencing of a minimal tiling path of cosmids, 
and, for selected areas, finished sequencing. During the 
construction of the transcript map, we identified new 
microsatellites that permitted further refinement of the 
candidate interval to the centromeric-most 115 kb. Us- 
ing two trapped exons as probes, we isolated a 3.7 kb 
clone from a leukocyte cDNA library. The protein prod- 
uct is predicted to be 781 amino acids in length and is 
a member of a family of proteins that includes the Ro52 
autoantigen (Itoh et al., 1 991 ) and RFP (ret finger protein, 
the gene encoded by the 5' half of the hybrid ret trans- 
forming gene; Takahashi et al., 1 988). We have identified 
three missense mutations in the cloned cDNA, all of 
which encode changes in the same highly conserved 
C-terminal domain. None of the three mutations is pres- 
ent in any of a large panel of normal control chromo- 
somes. Two of these mutations are found on FMF carrier 
chromosomes in populations that, in some cases, have 
been geographically separated for over 2000 years, indi- 
cating that most cases of FMF are descended from a 
very ancient pool of founders, and suggesting common 
origins for several Middle Eastern populations. 



Results 

Positional Cloning Strategy 

Figure 1 depicts schematically the recent steps in our 
positional cloning strategy. Based on early genetic local- 
ization data, we established a contig spanning the 1Mb 
interval between D16S94 and D16S2622 (Sood et al. t 
1997). Microsatellite markers from this interval defined 
a critical region of 285 kb (D16S468~D16$3376) based 
on observed recombinations in a panel of 61 families, 
and a smaller candidate interval of —200 kb {D16S3082- 
D16S3373) based on historical recombinants. Tran- 
scripts were then identified by a combined strategy of 
exon amplification, direct cDNA selection, and single- 
pass sequencing. This led to the isolation of 9 full-length 
cDNA clones (8 of which were previously uncharacter- 
ized), including 3 classical zinc finger genes, 2 members 
of the olfactory receptor family, and a novel caspase. 
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Figure 1. Schematic Representation of the MEFV Region on Chromosome I6p13.3 

(A) M Mb region over which a cosmid contig was constructed (Sood et ai., 1997). PKD1, polycystic kidney disease 1 ; TSC2, tuberous sclerosis 
2; RSTS, Rubinstein-Taybi syndrome; CREBBP, CREB binding protein gene. 

(8) --285 kb candidate interval defined by intrafamilial recombination events in our panel of families (Balow et af., 1997). D16S340S-D16S3373 
is the minima* candidate interval defined by historical recombinants in the present study. 

(C) Representation of the genomic structure of the v75-1 gene. Shaded boxes represent exonS; introns are drawn to scale. Numbers above 
the boxes represent the size of the exons, in bp. Numbers below the boxes reflect the order of the axons, with (1) being the most 5\ 



The furthest centromeric of these 9 cDNAs, clone v75-t , 
was isolated by solution hybridization of a leukocyte 
cONA library with biotinylated oligonucleotide probes 
derived from two exons trapped from PAC 273L24. 

In parallel, analyzing finished sequence from —90 kb 
in the center of the O16S3082-D16S3373 interval, we 
found two new microsatellites, D16S3404 and 0163340$ 
(Figure 1B). In one non-Ashkenazi Jewish family, evi- 
dence for a historical recombination event that occurred 
between O16S3404 and D16S3405 in the highly con- 
served non-Ashkenazi Jewish haplotype (designated 
haplotype A) was observed, thereby excluding the re- 
gion telomeric of D16S3405 (and four candidate genes 
encoded therein) from further consideration. 

Characterization of cDNA v75-1 

The translated v75-1 cDNA sequence is shown in Figure 
2, and the exon-intron structure deduced from the geno- 
mic sequence of two cosmids is depicted in Figure 1C. 
Although there is an excellent Kozak consensus (Kozak, 
1 996) at the initial methionine, the reading frame remains 
open in the cDNA upstream of this. Nevertheless, there 
are no splice-acceptor consensus sequences or in- 
frame methionines with good Kozak sequences before 
the first stop upstream in genomic DNA, and the tran* 
script size by Northern blot (see below) is 3.7 kb, thus 
suggesting that the sequence presented is full-length. 

Consistent with the initial analysis of single-pass ge- 
nomic sequence from this region, BLASTX (Altschul et 
al M 1990) demonstrated homologies with a number of 
molecules implicated in inflammation. These include the 
52 kDa Ro/SS-A ribonucleoprotein, against which pa- 
tients with systemic lupus erythematosus (SLE) and Sjo- 
gren's syndrome frequently make autoantibodies, the 



interferon-inducible transcriptional regulator Staf-SO 
(Tissot and Mechti, 1 995), and rpt-1 , a mouse down- 
regulator of lt-2 (Patarca et al., 1 988). Homology to the 
latter protein is in a domain extending from residues 
385 to 550, while homology to the former two proteins 
is particularly high through a domain at the C-termtnal 
end of v75-1 (double-boxed in Figure 2). This domain is 
also found in the ret finger protein (RFP) and has been 
termed an rfp (Takahashi et al., 1988) or 830.2 domain 
(Vemet et al., 1 993; Henry et al., 1997). Analysis with the 
SEG algorithm (Wootton, 1 994) indicated a high likeli- 
hood that this domain assumes a globular conformation. 

Additional analyses of the putative protein product, 
termed pyrin to connote its relationship to fever, indicate 
that it is very positively charged, with a predicted pi > 
8 and with lysine and arginine making up 13% of the 
amino acid composition. Moreover, the PSORT algo- 
rithm (Nakai and Kanehisa, T992) detected two overlap- 
ping nuclear targeting signals. The first is a four-residue 
pattern composed of a histidine and three lysines (start- 
ing at residue 41 9), and the second is a Robbins/Ding- 
wall consensus (boxed in Figure 2, residues 420-437; 
Robbins et al., 1991). In addition, a PROSITE search 
(Bairoch et al., 1 997) identified a bZIP transcription fac- 
tor basic domain (Shuman et al., 1 990) at residues 266- 
280. The region between residues 375 and 407 (denoted 
by plus signs in Figure 2) contains an arrangement of 
cystine and histidine residues with spacing that con- 
forms closely to a B box-type zinc finger domain (Reddy 
et al., 1992). 

Northern analysis using a probe with unique sequence 
from exon 2 of v75-1 showed an abundant ~>3.7 kb 
transcript in peripheral blood, but there was no signifi- 
cant expression in 14 other normal tissues, including 
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Figure 2. The Predicted Protein Sequence Encoded by cDNA v75-1 

The boxed segment from aa 266 to 2B0 is a bZIP transcription factor basic domain; the boxed segment from aa 420 to 437 is a flobbins/ 
Dingwalt consensus nuclear targeting signal. The segment indicated by plus signs between residues 375 and 407 is a potential B-box zinc 
finger domain. The region double-boxed from residue 577 to 757 is a rfp, or B30.2, domain. Within this box, the three mutated amino acids 
are double -underlined. 



spleen and thymus (Figure 3A}. Bone marrow and lymph 
nodes were also negative for v75-1 expression (data not 
shown). Among a panel of tumor celts, only the SW480 
colorectal adenocarcinoma cell line showed significant 
expression. Expression of v75-1 was absent in HL-60 
promyelocyte leukemia, K-562 erythroleukemia, and 
Burkitt's lymphoma. Even by RT-PCR, message was 
not detectable in a peritoneal fibroblast cell line, or in 
fetal or adult brain (data not shown). Figure 3B demon- 
strates pyrin expression in PMNs, but not in lympho- 
cytes. We also found that pyrin is not expressed in 
monocytes purified from peripheral blood (data not 
shown). 

Mutational Analyses 

Direct sequencing of RT-PCR products or amplified ex- 
ons from the eight cDNAs telomeric to v75-1 failed to 
identify disease- associated mutations. In contrast, three 
different missense mutations were identified in exon 1 0 
of v75-1 (Figure 4) after screening a total of 1 65 individu- 
als from 65 families. The first mutation, a G— C transver- 
sion at nt 2040 that results in the substitution of isoleu- 
cine for methionine (M680I), was observed in the 
homozygous state in the affected offspring of a single 
Armenian family. The second, an A^G transition at nt 



2080 causing a substitution of valine for methionine 
(M694V), was observed in a large number of affected 
individuals bearing four apparently distinct disease- 
associated haplotypes (see below). The third, a T— C 
transition at nt 21 77, results in the substitution of alanine 
for valine (V726A) and was observed in affected individu- 
als bearing the C haplotype in a Druze family and in 
other FMF patients and carriers bearing this haplotype. 
When the remaining nine exons of v75-1 were se- 
quenced in these haplotypes, no additional mutations 
were identified. 

Ancestral Relationships Among FMF 
Carrier Haplotypes 

Table 1 explores the relationship of v75-1 mutations with 
specific disease-associated haplotypes in our study 
population. Haplotype A is the disease-associated hap- 
lotype observed in a very high percentage of North Afri- 
can Jewish carrier chromosomes, but that is seen less 
frequently in the Armenian and Iraqi Jewish populations. 
Haplotypes F and G are additional North African Jewish 
disease-associated haplotypes, each observed in a sin- 
gle family, that appear by microsatellite analysis to be 
different from the A haplotype (Table 1A). Haplotype B 
is yet another distinct FMF haplotype observed in the 



Figure 3. Expression Profile of the v75-1 
Gene 
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noma SW480. HL-60, promyelocytic knjke- 
mia; K-56Z, erythroleukemia; MOLT4, iym- 
phoblastic leukemia; A549, lung carcinoma; 
G3S1 , melanoma. 

(B) Hybridization of the same exon 2 probe on 
Northern blot with mRNA from purified PMNs 
and lymphocytes. PMN lanes represent prep- 
arations from different individuals. A |i-actin 
control is shown below. 




Iraqi Jewish population, while the C haplotype is seen 
in patients from several different ethnic groups. The J 
and K haplotypes have only been observed in Armenian 
affected individuals. Table 1A indicates no easily dis- 
cernible relationship among these haplotypes at the 
level of microsatellites, except that haplotypes J and K 
are identical in the region from O16S3403 to D16S2617. 

Table 1B summarizes data on 11 single nucleotide 
polymorphisms (SNPs) that we have identified from 
within the v75-1 gene for each of the haplotypes defined 
in Table 1A. Nine of these polymorphisms represent 
synonymous substitutions within exons of this gene, 
while the remaining two are seen in introns 4 and 6. The 
A, F, and G haplotypes are identical at all 1 1 polymor- 
phisms, and, although the B haplotype differs from A, F, 
and G at the 5' end of the gene, SNPs for ail 4 haplotypes 
converge beginning at exon 3. Similarly, SNPs for the 
C, J, and K haplotypes also converge by exon 3. The 
"AFGB" and "CJK" core SNP haplotypes from exon 3 
to 1 0 were also observed on significant numbers of 
ethnically matched normal chromosomes (not shown). 

Table 1 C, which correlates haplotypes with v75-1 mu- 
tations, suggests that ancestral relationships among 
these subsets of chromosomes are likely. Thus, the con- 
vergent A, F, G, and B haplotypes all bear the M694V 
mutation, while we have not observed M694V on haplo- 
types that do not share the AFGB core SNPs. Similarly, 
Armenian haplotypes J and K converge within v75-1 and 
both bear the same M680I mutation, while this mutation 
has not been observed on haplotypes that do not have 
the CJK core SNPs. The C haplotype bears the same 
convergent SNP haplotype as J and K, suggesting an 
independent mutational event on an ancestral CJK-type 
chromosome. 



For any given haplotype, only one mutation was pres- 
ent. Altogether, these three mutations accounted for 78 
carrier chromosomes in our panel. There were several 
other carrier chromosomes with different microsatellite 
haplotypes for which we did not observe any of these 
three mutations (these latter carrier chromosomes have 
thus far only been screened for exon 10 mutations). In 
four separate families we have observed M694V/V726A 
compound hetero zygotes. In one such family, the par- 
ents are of Ashkenazi and Iraqi Jewish ancestry, with 
the Ashkenazi father carrying the V726A mutation on 
the C haplotype FMF chromosome. This haplotype es- 
tablishes an ancestral link between Ashkenazi Jews, 
who have lived in Eastern Europe for the last 1 500 years, 
and the Druze, who have lived in isolated communities 
in what is now Israel and Lebanon for the last —1000 
years. The only other Ashkenazi Jewish carrier whom 
we have genotyped also has the V726A mutation. 

The M680I, M694V, and V726A mutations were not 
present on any of 1 05 normal chromosomes in our study 
population (defined by haplotype analysis in well -char- 
acterized families), nor in almost 200 North American 
control chromosomes. It is also noteworthy that al- 
though no A haplotype— bearing noncarrier chromo- 
somes could be identified, two noncarrier chromosomes 
with the B haplotype (both by microsatellites and SNPs) 
did not have the M694V mutation, further strengthening 
the connection between the presence of the mutation 
and disease. Moreover, 42% of a panel of 24 Middle 
Eastern normal chromosomes bear the AFGB SNP core 
haplotype but do not carry M694V, which is further evi- 
dence that M694V is disease-causing, and not a private 
polymorphism. 
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Figure 4. DNA Sequence Electrapherograms Demonstrating the 
M680I, M694V, and V726A Substitutions 

Fof each mutation, individuals who are homozygous for the normal 
allele are shown at the top, heterozygotes between the normal and 
mutant allele are shown in the middle, and homozygotes for the 
mutation are shown at the bottom. On the pedigree diagram, the 
heterozygote is shown as an open red symbol, and the hornozygote 
for the mutation is shown as a closed red symbol. The control for 
the M680I mutation is individual 76-06, the control for M694V is 76- 
02, and the control for V726A is 3*01 . 



Discussion 

In this paper, we describe the molecular cloning of a 
new cDNA, v75-1 , from the familial Mediterranean fever 
candidate region on chromosome 16p13.3 and present 
compelling evidence that this gene is MEFV, the gene 
causing FMF. The cornerstone of the argument equating 
v7S-1 with MEFV is the identification of three different 
v75-1 mutations on FMF carrier chromosomes in multi- 
ple ethnic groups, mutations that are not seen in a panel 
of almost 300 normal control chromosomes. 

It is extremely unlikely that the substitutions we have 
identified in v75-1 are actually polymorphisms in tight 
linkage disequilibrium with "real" mutations on a nearby 
gene. This hypothesis would require that there be three 
such v75-1 polymorphisms on three different haplo- 
types, each in perfect linkage disequilibrium with the 
mutations on the real FMF gene. While not impossible, 
such a scenario is at the least unnecessarily complex. 
It is also unclear where such a closely linked gene would 
be located. The historical recombinants at the 5' (centro- 
meric) end of v75-1 exclude the interval between 
D16S3373 and v75-1. On the telomeric side, the 5' end 
of a novel zinc finger gene is located within 1 0 kb of the 
3' end of v75-1 , but thorough screening has revealed no 
mutations in this latter gene (data not shown). Moreover, 
there are no trapped exons, direct-selected cONAs, or 
expressed sequence tag (EST) hits that map to the inter- 
val between them. Finally, and most importantly, the 



observation of normal chromosomes that bear disease- 
associated microsateilite and SNP haplotypes but do 
not have the M680I, M694V, or V726A mutations is strong 
evidence that these are not just haplo type- specific poly- 
morphisms. 

While these genetic arguments firmly establish the 
identity of v75-1 as the FMF gene, data on v75-1 tissue 
distribution are highly consistent with the clinical pheno* 
type. Based on the nature of the inflammatory infiltrate 
and the anatomic localization of inflammation in FMF, 
MEFV gene expression might be predicted to be ob- 
served in granulocytes and/or serosal cells. Multiple- 
tissue Northern blots demonstrated high levels of ex- 
pression in peripheral blood leukocytes, but not in lymph 
node, spleen, or thymus, which are comprised largely 
of lymphocytes. From this we deduced that pyrin is most 
likely expressed primarily in mature granulocytes, an 
inference that was subsequently confirmed by direct 
observation. 

Computational analyses of wild-type and mutant pyrin 
molecules suggest a mechanism by which the mutations 
identified may cause disease. AN three mutations are 
clustered within 46 amino acids of one another in the 
highly conserved B30.2 (rfp) globular domain at the 
C-termina! end of the predicted protein, and they proba- 
bly affect the secondary structure of this domain. Based 
on BLASTX homologies and its restriction to granulo- 
cytes, pyrin is likely to be a nuclear factor that controls 
the inflammatory response in differentiated PMNs. 
Given that FMF is a disease of excessive inflammation, 
and that pyrin is homologous to rpt-1 , a known down- 
regulator of inflammation, pyrin may be a negative auto- 
regulatory molecule in PMNs. The FMF-associated mu- 
tations and resultant pyrin structural changes would 
then, after some undefined event that triggers inflamma- 
tion, prevent the norma! pyrin -mediated negative feed- 
back loop. 

None of these three mutations results in a truncated 
protein, and the periodic nature of inflammatory attacks 
in FMF is consistent with a protein that functions ade- 
quately at steady state but decompensates under 
stress. Several other diseases with periodic manifesta- 
tions follow this paradigm, including sickle cell anemia 
(p-gtobin; Weatherall et al., 1 995) and hyperkalemic peri- 
odic paralysis (adult skeletal muscle sodium channel; 
Ptacek et al., 1 991 ). Both are characterized by missense 
mutations encoding an abnormal protein with reduced 
function. Patients experience symptoms when environ- 
mental factors (hypoxia or acidosis in the case of sickle 
cell, excess consumption of potassium in periodic paral- 
ysis) induce additional conformational changes in the 
protein (sickle cell) or destabilize a delicate equilibrium 
(periodic paralysis). Similarly, environmental triggers 
may play a role in FMF, and the identification of pyrin 
should provide clues regarding the proximal steps in 
the cascade that leads to an attack. 

Just as protein-truncating mutations in globin pro- 
duce thalassemia and not sickle cell disease, more pro- 
found loss-of -function mutations in pyrin may produce 
a phenotype much different from FMF. It is even conceiv- 
able that the total lack of pyrin may lead to generalized, 
uncontrolled inflammation and may therefore be embry- 
onically lethal. Mouse models may shed light on these 
questions. 
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Table 1 . FMF Carrier Haplotypes and cONA v75-1 Mutations 
A. Extended Microsatellite Haplotypes 



Microsatellite Markers' 



Hap)otype 


Ethnic Grp° 


S468 


S3062 


$3403 


S3405 


S3404 


S3370 


S26T7 I 


S3373 


S3275 


S3376 


A 


NAfJ, IrJ, Arm 


212 


148 


1100 


165 


118, 122 


153 


104 


185 


184 


211 


F 


NAfJ 


204 


152 


1000 


165 


118 


149 


104 


187 


188 


215 


G 


NAfJ 


202 


150 


1000 


185 


94 


141 


107 


187 


184 


211 


B 


IrJ 


206 


138 


1000 


187 


94 


141 


110 


185 


186 


217 


C 


IrJ, Ash J, Dr, Arm 


202 


140 


1100 


165 


110 


149 


107n c 


187 


186 


201 


J 


Arm 


202 


136 


1100 


165 


118 


149 


107n 


183 


188 


215 


K 


Arm 


202 


132 


1100 


165 


118 


149 


107n 


185 


196 


211 



B. cONA v75-1 Intragenic Haplotypes 

Single Nucleotide Polymorphisms" 







x2.1 


■< 2.2 


<2.3 


x2.4 


x3.1 




x5.1 


:<5.2 


*5.3 




x9.1 


Haplotype 


Ethnic Grp 


306 


414 


495 


605 


901 


i4.1* 


1422 


1428 


1530 


16. 1* 


1764 


A 


NAJ, IrJ, Arm 


C 


G 


A 


A 


C 


A 


G 


A 


T 


c 


A 


F 


NAfJ 


C 


G 


A 


A 


c 


A 


G 


A 


T 


c 


A 


G 


NAfJ 


c 


G 


A 


A 


c 


A 


G 


A 


T 


c 


A 


B 


IrJ 


T 


A 


C 


G 


c 


A 


G 


A 


T 


c 


A 


C 


trj, AsrW, Or, Arm 


T 


A 


C 


G 


T 


G 


A 


G 


C 


T 


G 


J 


Arm 


T 


A 


C 


G 


T 


G 


A 


G 


C 


T 


G 


K 


Arm 


C 


G 


A 


G 


T 


G 


A 


G 


c 


T 


G 



C. cDNA v75-1 Mutations 



Mutation aa Haplotypes Norma) Chromosomes 



2040 G—C Met680lle J (1/1) 


0/293 


K(1/1) 




2080 A— G Met694Val A (62/62) 


0/295 


B (5/5) 




F(1/1) 




G(1/1) 




21 77 T— C Val726Ala C (8/8) 


0/283 


•In each case the "0 16" has been omitted at the beginning of the locus designation. The position of MEFV relative to these markers is indicated 


by the arrow. Microsatellrtes are listed tetomeric to centromeric, from left to right. 


"NAfJ = North African Jewish, AshJ = Ashkenazi Jewish. IrJ = iraqi Jewish, Arm = Armenian, Or = Oruze. 


'Denotes the presence of a null allele with standard D16S261 7 primers, but 107 bp amplicon when an alternative primer pair is used (Salow 


et al., in press). 




"Designated by exon or intron; <2.1 indicates the first polymorphism in axon 2. Below each designation is the nucleotide position in the 


sequence, with 1 being the adenosine in the ATG initiation codon. Note that the polymorphisms are listed 5' to 3', from left to right, which is 


the opposite orientation on the chromosome from that shown in (A). 




♦Intron 4, nucleotide position 1356 - 43. 




'Intron 6, nucleotide position 1610 + 95. 





It would appear that at least some of the phenotypic 
variation in FMF may be attributable to the difference 
between the M694V and V726A mutations. The M694V 
mutation is very common in the populations with the 
highest incidence of systemic amyloidosis (especially 
North African Jews), while V726A is seen in populations 
in which amyloid is less common (Iraqi and Ashkenazi 
Jews, Druze, and Armenians). Nevertheless, M694V- 
associated haplotypes can comprise a significant mi- 
nority of carrier chromosomes in populations where se- 
vere amyloidosis is exceedingly rare. One possible 
explanation could be that the "milder" V726A mutation 
is "protective" against amyloid, in which case M694V7 
M694V hompzygotes would be amyloid-prone but 
V726A homozygotes and compound heterozygotes 
would not. 

Although founder effects clearly have played a role in 
the population genetics of FMF, selection for heterozy- 
gotes remains an attractive hypothesis. The probable 



rote of pyrin in the homeostatic control of inflammation 
is consistent with heterozygotes exhibiting heightened 
responses to a specific pathogen or class of pathogens. 
The strikingly high gene frequencies for FMF in more 
than one Middle Eastern population could then be ex- 
plained by increased heterozygote resistance to such 
an organism. Since the gene frequency for FMF is much 
lower in the East European Ashkenazi Jewish population 
(Yuval et al., 1 995), it is possible that such a selective 
factor would have been endemic to the Mediterranean 
basin. The identification of pyrin should lead to testable 
hypotheses regarding the host inflammatory response 
in FMF heterozygotes. 

Our data indicate that the M694V and V726A muta- 
tions are very old. Based on the convergence of several 
carrier haplotypes, it is likely that all modem M694V 
chromosomes are descended from a common founder. 
This mutation is seen on the B FMF haplotype commonly 
observed in Iraqi Jews, a group that remained relatively 
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expressed sequence tag (EST) hits that map to the inter- 
val between them. Finally, and most importantly, the 



onically lethal. Mouse models may shed light on these 
questions. 
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Figure 5. Multiple Sequence Alignment of 
the C-Terminal End of Human Pyrin 
Sequences shown have statistical similarity 
as assessed by BLAST (Aitschul et al., 1 990). 
Search cutoffs used to identify homologs 
were a Kartin-Altschul score of two aligned 
sequences 70 with a probability ^10 '.At 
each position, residues occurring in a major- 
ity of the sequences are shown tn inverse 
type. The numbering scheme at the top of the 
figure is based on the sequence of pyrin. The 
positions of secondary structural elements in 
the wild-type pyrin (here, all 3 sheets) as pre- 
dicted by the profile neural network method 
PHOsec (Rost and Sander, 1993, 1994) are 
shown as bold, horizontal arrows below the 
alignment. The residues altered in the three 
mutants are indicated by downward -pointed 
arrowheads above the alignment. hum-RFP. 
RET finger protein {SWISS-PROT P14373); 
xia-xnf7, nuclear phosphoprotein xnf7, Xeno- 
pus laevis (PIR A43906); pwa-A33, zinc-bind- 
ing protein A33, Pteurodeles waltl (SWISS- 
PROT Q02084); hum-SS-A/Ro, 52 kDa RO 
protein (SWISS-PROT P19474); hum-afp, acid 
finger protein (GenSank U09825); hum-BT, 
butyrophilin (GenSank U90552); hunvefp, es- 
trogen-responsive finger protein (PIR A49666); 
hum-B30-2. B30-2 gene, (PRF 2002339); pig- 
RFB30, ring finger protein RFB30. Sus scrofa 
(EMBL Z97403); hum-Staf-50, transcription 
regulator Staf-50 (PIR A57041). 
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isolated from the rest of the Jewish population from the 
time of the Babylonian captivity (-2500 years ago) until 
the formation of the state of Israel. Similarly, the obser- 
vation of the V726A mutation and C microsatellite haplo- 
type in Armenians, Ashkenazi Jews, and Druze FMF 
patients establishes a historical lower boundary on the 
age of this mutation of perhaps 2000 years. 

It is also possible, of course, that the M694 V mutation, 
and its surrounding SNPs, has been transferred to four 
different hapiotypes by interchromosomal gene conver- 
sion rather than by meiotic recombination events. In that 
case, the estimation of the age of the original founder 
mutation would not be possible. Finally, the possibility 
of multiple independent origins of M694V and M680I 
cannot be formally excluded, but the concordance of 
the SNP hapiotypes close to the mutation (Table 1 B and 
1 C) renders this possibility rather unlikely. 

The identification of pyrin as the molecule responsible 
for FMF also focuses interest on the RoRet family of 
proteins, all of which share a B-box zinc finger as well 
as a 170 amino acid rfp (830,2) domain. The spacing 



between these two structural elements is also highly 
conserved, suggesting that precise orientation of the 
two domains with respect to one another is required for 
function. The conservation of the rfp domain (Figure 5) 
in molecules as diverse as butyrophilin (a milk protein 
with probable receptor function; Jack and Mather, 1 990), 
Ro52 (a ribonucleoprotein that is autoantigenic in pa- 
tients with SLE and Sjogren's syndrome), A33 (a factor 
that binds polytene chromosomes in the newt; Bellini 
et al., 1 993), and xnf7 (a factor that binds mitotic chromo- 
somes in the frog; Reddy et aJ., 1991) is remarkable; the 
sequence identity among these proteins ranges from 
40% to 60%. Scrutiny of B30.2 mutations in pyrin may 
provide insights into this evolutionariiy conserved but 
poorly understood domain. 

The identification of pyrin mutations as the cause of 
FMF will substantially advance our understanding of this 
disorder and will provide a tool to address questions 
ranging from population genetics to structural biology. 
It represents the first step in the delineation of an impor- 
tant new pathway in the control of inflammation. 
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Experimental Procedures 
Patients and DNA Samples 

Forty-four families of non-Ashkenazi Jewish descent {18 Moroccan, 
14 Libyan, 5 Tunisian. 2 Egyptian, and 5 Iraqi) and 5 Arab/Druze 
families were identified and sampled at the Chaim Sheba Medical 
Center in Tef-Hashomer, Israel. Twelve Armenian families were re- 
cruited from Cedars-Sinai Medical Center in Los Angeles. Details 
of this panel have been previously described (Balow et al., 1997). 
One additional Ashkenazi/iraqi Jewish family was also studied. The 
diagnosis of FMF in all families was according to established clinical 
criteria (Soharet al.. 1967). Patient consent was obtained from parti- 
cipants in the study after approval by the human experimentation 
committees at each institution. DNA was extracted from whole blood 
or from Epstein- Barr virus-transformed lymphocytes by standard 
techniques. 

Genotyping 

Patients were genotyped as previously described (Balow et al., 
1997). Additional, genotyping by PCB was carried out with three 
new markers that were developed in this study and are described 
here (D16S3403, D16S3404, and 016S3405). Sequences of oligonu- 
cleotides and PCFt conditions used to amplify these markers are 
available from the Genome Database. 

Direct cONA Selection 

Cosmids, BAC, and P1 clones in the FMF candidate region were 
biotinytated using BioPrime (Life Technologies, Gaithersburg, MD). 
cONAs were prepared from combined mRNA from fetal brain, fetal 
liver, and human lymph node by reverse transcription and ligation 
of an EcoRI/Notl adaptor, which also served as a PCR primer, to 
second-strand cONAs. cDNAs were directly hybridized to biotinyl- 
ated templates that were recovered using streptavidin-labeled mag- 
netic beads. Conditions for blocking, hybridization, binding, and 
elution of cDNAs from magnetic beads (Dynal) were as described 
(Parimoo et al., 1991). After two rounds of selection, eluted cDNAs 
were amplified with CUA-tailed EcoRI/Notl adaptor primers and sub- 
cfoned into the pAMPlO vector (Life Technologies, Gaithersburg, 
MD) to yield libraries of selected cDNAs. Recombinant clones were 
arrayed on blots. Clones that hybridized to either repetitive or ribo- 
somai sequences were excluded from further analysis. To confirm 
their origin, unique clones were individually hybridized to EcoRI 
digests of cosmid/BAC/P1 DNAs and DNAs from chromosome 1 6- 
specific human-hamster hybrid lines. Clones were then hybridized 
to each other and were binned into groups. Representative clones 
of each group were hybridized to multiple-tissue Northern blots and 
sequenced, 

Exon Trapping 

Exon trapping was performed on PAC clone 273L24 (Sood et al,, 
1997) as previously described (Buckler et al., 1991). PAC clone 
273L24 was partially digested with Sau 3A1. The reaction products 
were size-fractionated by agarose gel electrophoresis, and DNA 
fragments 2 kb and larger were isolated from the gel. Partially di- 
gested DNA (50 ng) was ligated with 10 ng of exon-trapping vector 
pSPL3 that had been previously cleaved with BamHI and dephos- 
phorylated with calf intestinal alkaline phosphatase (Promega, Madi- 
son, Wl). Ligation products were electroporated into E. coli OH1 2B, 
and the electroporated cells cultured en mass in LB broth with 200 
mg/ml ampicillin for 1 6 hours at 37 C C with shaking. ONA prepared 
from the culture was used to transfect COS-7 cells (ATCC 30-2002) 
using lipofectACE reagent (Life Technologies, Gaithersburg, MD). 
Total RNA was isolated from transfected COS-7 cells with Trizot 
reagent (Life Technologies) followed by ethanol precipitation. First- 
strand cONAs of transcription products from pSPL3 were primed 
with the oligonucleotide SA2. Specific amplification of trapped ax- 
ons was as follows: PCR primed with oligonucleotides SA2 and SD6 
was performed, followed by digestion of the PCR products with 
BstXt. A second PCR reaction using the digestion products was 
primed with oligonucleotides dUSD2 and dUSA4. The resulting DNA 
fragments were cloned into pAMP10 vector and sequenced. Two 
hundred clones were sequenced and 20 independent exons were 



identified by visual inspection and hybridization to DNA fragments 
from the FMF critical region, with several exons identified more than 
one time. 

cDNA Identification by Solution Hybridization 
Solution hybridizations were carried out using the GeneTrapper 
cDNA Positive Selection System (Life Technologies, Gaithersburg, 
MD). Two trapped exons, v66 and v75, were used as starting mate- 
rial. PCR screening of Superscript c DNA libraries (Life Technologies, 
Gaithersburg, MD) derived from human brain, liver, leukocytes, 
spleen, and testis were used to determine the tissue-specific ex- 
pression of these exons. The leukocyte cDNA library was used in 
further experiments since it showed significant expression of v66 
and v75 exons. GeneTrapper experiments were performed with 
sense and antisense primers from both exons, assuming both orien- 
tations of these exons in the putative transcript. The following oligo- 
nucleotides were synthesized and PAGE- purified: v66GTl : AAG CTC 
ACT GCC TTC TCC TC: V66GT2: GAG GAG AAG GCA GTG AGC 
TT; V7SGT1 : GAC TTG GAA ACA AGT GGG AG; V75GT2: CTC CCA 
CTT GTT TCC AAG TC. Oligos were biotinylated and hybridized to 
single-stranded DNA from the leukocye cONA library (one primer 
per reaction), followed bycDNA capture using paramagnetic strep- 
tavidin beads and repair using the corresponding nonbiotinylated 
oligos. Colony hybridization of lifts using n P-dCTP end-labeled 
oligos was used to identify positive clones- Gel-purified inserts from 
these clones were hybridized to cosmid contig blots to distinguish 
cDNA clones mapping to the FMF region from false positive clones 
due to homologous domains. All positive clones were identified by 
the primers V66GT2 and v76GT2, and no clones were identified by 
the other set of primers. 

Oligonucleotides for Exon Amplification 

Oligonucleotides used to amplify pyrin exons were as follows (all 
oligo sequences are given 5' to 3'}: exon 1 forward, AACCTGCCTTTT 
CTTGCTCA; exon 1 reverse, CACTCAGCACTGGATGAGGA; exon 3 
forward, GAACTCGCACATCTCAGGC; exon 3 reverse, AAGGCCCA 
GTGTGTCCAAGTGC; exon 4 forward, TTGGC AC CAGCTAAAGAT 
GGC; exon 4 reverse, TCTCCCTCTACAGGGATGAGC; exon 5 for- 
ward, TATCGCCTCCTGCTCTGGAATC; exon 5 reverse, CACTGTGG 
GTCACCAAGACCAAG; exon 6 forward, TCCAGGAGCCCAGAAGTA 
GAG; exon 6 reverse. TTC TCC C TATC AAATCCAG AG ; exon 7 
forward, AG AATGTAGTTC A TTTCC AGC; exon 7 reverse. CATTTCTG 
AACGCAGGGTTT; exon 8/9 forward, ACCTAACTCCAGCTTCTCTC 
TGC; exon 8/9 reverse, AGTTCTTCTGGAACGTGGTAG: exon 10a 
forward, CCAGAAGAACTACCCTGTCCC; exon 10a reverse, AG AGC 
AGCTGGCGAATGTAT-, exon 10b forward, GAGGTGGAGGTTGGAG 
ACAA; exon 10b reverse. TCCTCCTCTGAAATCCATGG. 

Northern Blot Analysis 

To determine transcript size and level of expression in various tis- 
sues, multiple-tissue Northern blots (Clontech) were hybridized with 
probes derived from various exons of the gene. These exons were 
amplified and purified as part of the sequencing protocol for muta- 
tion analysis. Larger exons (2, 5, and 10) were labeled by random- 
priming using Stratagene Prime-It Kit and )2 P-dCTP (ICN). Hybridiza- 
tion and washing of blots were essentially as described in Sambrook 
et al. (1 989}, except using Hybridisol I (Oncor) prepared hybridization 
buffer. Hybridization was detected by autoradiography, with 4 hour 
exposures. Northern blots with mRNA from highly purified peripheral 
blood lymphocytes, PMNs, and monocytes were the kind gift of Ore, 
H. Lee Tiffany and Harry Malech. 

Mutation Detection by Fluorescent Sequencing 
Approximately 1 00 ng of genomic DNA template was used in PCR 
reactions to amplify exons and flanking intronic sequences ac- 
cording to the supplier's recommendations for AmpliTaq Gold (Per- 
kin Elmer, Branchburg, NJ) and Advantage-GC Genomic PCR Kit 
(Clontech, Palo Alto, CA). The PCR primers were tailed with one of 
the following sequences: -21M13 forward: GTAAAACGACGGC 
CAGT; -28 M1 3 reverse: CAGGAAACAGCTATGACCAT; -40 M1 3 for- 
ward*. GTTTTCCCAGTCACGACG, After amplification, reactions 
were run on t% agarose gels and gel-purified using either QiAquick 
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gel extraction kit (QIAGEN, Santa Clarita, CA) or Microcon/Micro- 
pure/Gel Nebulizer system (Amicon, Beverly, MA). Alternatively, PCR 
products were column-purified with Microcon-100 (Arnicon). Puri- 
fied amplicons were sequenced with dye primer chemistry (PE Ap- 
plied Biosystems, or Amersham, Cleveland, OH). Sequencing reac- 
tions were ethanot -precipitated and njn on an AB1 377 automated 
sequencer. Sequence data were analyzed with either Autoassem bier 
1 .4 (PE Applied Biosystems, Branchburg, NJ) or Sequencher 3.0 
(Gene Codes inc., Ann Arbor, Ml). 
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